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ABSTRACT 

We present near-infrared spectra obtained to search for Paa and molecular hydrogen 
lines in edge-darkened (FR I-type) radio galaxies with bright Ha emission in the redshift 
range 0.0535 < z < 0.15. We find that all three galaxies in our sample (PKS 0745-191, 
PKS 1346+26, & PKS2322-12) which are associated with strong cooling flows also have 
strong Paa and H2l-0S(1) through S(5) emission, while other radio galaxies do not. 
Together with earlier observations this confirms claims that cooling flow galaxies are 
copious emitters of molecular hydrogen with large H2 1-0S(3)/Paa ratios in the range 
0.5 to 2. The emission is centrally concentrated within the inner few kiloparsec and could 
come from warm (T~1000-1500 K) molecular material which is being deposited by the 
cooling flow. We speculate that the H2 emission could be related to the interaction 
between the jets and this molecular gas. 

Subject headings: galaxies: clusters: general - cooling flows - galaxies: ISM - galaxies: 
jets - galaxies: ISM - galaxies: elliptical and lenticular, cD 
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1. Introduction 



2. Observations 



Cooling flows (hereafter CFs) arc usually invoked 
to explain the huge, extended X-ray himinosity seen 
in many massive elliptical galaxies and clusters. They 
are thought to represent a massive, more or less spher- 
ical, subsonic inflow of hot gas due to radiative cool- 
ing in the central regions of the galaxy (e.g. Fabian 
1994). The typical inflow rates are of the order of 
several tens to several hundreds of solar masses per 
year. However, it has been a puzzle where the huge 
amount of inflowing matter is actually deposited after 
it has cooled, since there was little direct evidence in 
any other waveband for this material (e.g. CO or H I, 
O'Dea et al. 1994, 1995). 

Elston & Maloney (1994) presented a brief report 
that strong 1-0 S(l) A2.1213/Ltm molecular hydrogen 
emission occurs in CF galaxies, however no individual 
detections are presented in that paper. Jaffe & Bre- 
mer (1997, hereafter JB97) also have detected strong 
H2 1-0 S(l) emission in three CF galaxies (Abell 478, 
Hydra-A, PKS0745-19) and none in a control sample 
of five galaxies, suggesting that such emission may be 
characteristic of CF galaxies. JB97 suggested that 
the molecular hydrogen emission lines, which usually 
probe shocked or X-ray heated gas with a temperature 
of 1000-2000 K, trace some of this missing matter. 

In this letter we report near-infrared spectroscopy 
of more CF and control galaxies which for the first 
time also include transitions from S(2) to S(5). The 
galaxies were observed in a program to search for hid- 
den, broad Paa lines in FRI radio galaxies, to test 
the idea that the optical and radio properties of FR I 
and FR II radio galaxies differ not because of a differ- 
ent central engine, but because of a 'closed' obscuring 
torus in FRI's (Falcke, Gopal-Krishna, & Biermann 
1994). For this program, we selected a number of 
bright radio galaxies with FRI or amorphous radio 
structure and luminous Ha emission in the redshift 
range 0.0535 < z < 0.15, so Paa would be redshifted 
into the K-band window. The galaxies were selected 
from the samples of Morganti, Ulrich, & Tadhunter 
(1992) and Owen et al. (1995). Since CF galaxies 
tend to be overluminous in emission lines compared 
to normal radio galaxies (e.g. Baum 1992), our sam- 
ple had a number of such objects. However, since the 
sample was selected with another objective, it also 
includes suitable control galaxies with similar proper- 
ties to those with CFs. 



Our spectra were obtained during two observing 
runs in April and December, 1996, at the Multiple- 
Mirror-Telescope (MMT) (operated jointly by the 
Smithsonian Astrophysical Observatory and the Stew- 
ard Observatory) with the FSPEC near-infrared long- 
slit spectrometer (Williams et al. 1993). Most of the 
time the weather was not photometric. The galaxies 
in our sample were observed in the K-band with either 
a 75 1/mm (2nd order) or a 300 1/mm (1st order) grat- 
ing. The slit-width was ~ 2.5 pixels (0'.'43pixel~-^), 
and guiding utilized an infrared camera viewing the 
spectrometer slit. The low-resolution grating covered 
the whole K-band range, while the wavelength range 
for the high-resolution grating was chosen such that 
it covered the range between Paa and H2l-0S(2). 
PKS 0745—191 was observed at two different central 
wavelengths to obtain a broader spectral coverage. 

During the observation, sets of four exposures of 
the galaxies were taken, each of which was offset by 
^7" along the slit. Every 60 90 mins we observed 
a nearby calibration star in the same manner. The 
spectra were sampled on the chip by 2-3 pixels per 
resolution element along the spatial and the disper- 
sion axes. The loci of the spatial and dispersion axes 
in each dataset were determined by tracing the OH 
airglow lines and the continuum of the calibration star 
in each of its four positions along the slit. The expo- 
sures were then geometrically transformed so that the 
spatial and dispersion axes became parallel with the 
axes of the array. 

The data were first reduced following standard pro- 
cedures for FSPEC (e.g. Engelbracht et al. 1996). We 
then re-reduced the data with a set of newly written 
IRAF routines which allow an almost fully automated 
data reduction and guarantee that all data sets, and 
specifically the OH airglow subtraction, are handled 
in an identical manner. The general approach was 
first to subtract the combined preceding and subse- 
quent frames from the raw-image, which effectively 
removes the dark current and most of the OH airglow 
lines. Each background subtracted exposure was then 
flatfielded, median smoothed with a box size com- 
parable to the instrumental resolution to correct for 
cosmic rays, and geometrically corrected. We then 
centered the continuum in a given aperture and ex- 
tracted a spectrum of the object and of the residual 
OH airglow ^14" offset from the galaxy in the same 
frame. Finally, the residual-subtracted spectra were 
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co-added with appropriate scales and weights, divided 
by the caUbration star (normaUzed to a black body 
spectrum) and wavelength calibrated using the OH 
airglow spectrum. We found that the automated pro- 
cedure yields spectra similar to those from the stan- 
dard procedures. Comparison of the two independent 
reductions gives assurance that the spectra are free of 
severe artifacts. 

3. Results 

The resulting spectra of three of our seven galaxies 
are shown in Figure 1. Due to the weather condi- 
tions during our observations and the resulting lack 
of good photometric data, we did not attempt a flux 
calibration and the flux densities are given in detec- 
tor units. For each spectrum we fitted the contin- 
uum with a polynomial and added/subtracted the one 
sigma error per pixel derived from the variations seen 
in the individual exposures during the co-adding pro- 
cess (dotted lines in Fig. 1). This procedure gives a 
robust estimate of the noise level in our spectra which, 
due to atmospheric absorption troughs and strong OH 
airglow lines, can be a function of wavelength. The 
pixel-by-pixel error estimates allow one to distinguish 
artifacts due to improper sky-subtraction from real 
spectral features. Parameters of the line identifica- 
tions are given in Table 1. Where line strengths are 
poorly measured, either due to terrestrial atmospheric 
interference or because the line fell near the end of the 
spectrum, the results are placed in parentheses in the 
table. 

In three cases (PKS 0745-191, PKS 1346-1-26, & 

PKS 2322-12), all strong CF galaxies, we have dis- 
covered strong Paa and molecular hydrogen lines. We 
did not find any evidence for either Paa or H2 in 
Abell 610, B2 0836-^29, or B2 0915+32. Although 
these galaxies have relatively luminous Ha emission, 
they are not known to have strong CFs. 

The most extreme emitter of molecular hydrogen 
is PKS 2322-12, where the H2 1-0 S(3) line has twice 
the flux of Paa and we also detected H2l-0S(4). 
In the other two CF galaxies, PKS 0745-191 and 
PKS 1346-^26, the H2 1-0S(3)/Paa ratio is roughly 
one half. PKS 0745—191 also shows emission from the 
H2 10 S(5) line at the end of the spectrum. With the 
exception of this line, whose flux is very uncertain due 
to poor determination of the continuum, all the line 
ratios we observe are consistent with thermal excita- 
tion being the dominant process, with a temperature 



T = 1000 1500 K. This temperature is similar to that 
found by Mouri (1994) for NGC 1275, another CF 
galaxy. 

It is known that CFs result in a detectable level of 
Ha emission, thought to be excited by shocks or pho- 
toionization (Heckman et al. 1989). However, since 
the differential extinction between Paa and the H2 
lines is negligible, this ratio is a more reliable indica- 
tor of excitation and the relative strength of molecular 
hydrogen emission than H2/Ha. Nevertheless, here 
we estimate the H2 1-0 S(l)/Ha ratios for our objects 
in order to facilitate a comparison with JB97. Assum- 
ing Ha/Paa = 5 (Hill, Goodrich & DePoy 1996) and 
a typical ratio 1-0 S(3)/l-0 S(l) = 0.5, we determine 
H2 1-0 S(l)/Ha « 0.2 (up to 0.8 for the extreme case), 
in agreement with JB97. 

From the width of the OH airglow lines, we deter- 
mined that our final spectral resolution (FWHM) in 
the 300 1/mm grating spectra is ~ 30 A (410kms~^ 
at 2.2/im). The line widths we mcasiire are clearly 
broadened with an observed FWHM of 40-50 A, which 
implies an intrinsic line width of several hundred 
kms~^ (after subtraction of the instrumental reso- 
lution in quadrature). 

For the 75 1/mm grating we achieved a resolution 
(FWHM) of - 60 A (820kms-i at 2.2/^m) and the 
lines in PKS 1346-1-26 are only barely resolved. Only 
the width of Paa seems to be large; however, this 
is probably an artifact due to the combination of a 
strong OH line and a deep atmospheric absorption 
trough around 2.0/Lim as indicated by a sudden in- 
crease of the noise in the red wing of the line. Also in 
PKS 0745—191 we find a suspicious, extremely broad 
base to the Paa line; however, there is a very good 
chance that this is due to a conspiracy of noise peaks. 

The lines in PKS 2322-12 and PKS 0745-191 
seem to be concentrated in the inner 2-3" of the 
galaxy (FWHM), while the emission in PKS 1346-f26 
seems to be more extended (~ 4-5"). Due to the 
limited signal-to-noise ratio of the spectrum of the 
low surface brightness emission in this galaxy and the 
small stepping increments, which lead to a spatial cut- 
off, it is difficult to specify the extent more precisely. 
We also do not see a variation of the H2 1-0S(3)/Paa 
ratio within the inner 3-4" of our galaxies exceeding 
15 per cent, which may indicate that the lines are 
co-spatial. 
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4. Summary &: Discussion 

Our detection of molecular hydrogen in two new 
CF galaxies and of additional strong H2 lines in a 
third, confirm the claims by Elston & Maloney (1994) 
and JB97 that CF galaxies are copious emitters of 
near-infrared molecular hydrogen lines. Moreover, 
this claim is further strengthened by additions to the 
list of control galaxies without CPs where H2 was not 
detected. 

There are now seven CF galaxies (NGC 1275, 
Abell 478, PKS 0745-191, Hydra A, PKS 1346+26, 
Cygnus A, & PKS 2322—12) which share the same 
properties: they have a) high estimated inflow rates 
of several lOOMQyr"^, b) relatively strong radio 
sources, c) luminous Ha emission, and d) luminous 
H2 emission (Table 2). The last property makes these 
galaxies outstanding — no other class of radio galax- 
ies has been found to have such a strong molecular 
hydrogen emission. In particular, of a total of eight 
(JB97 and this paper) galaxies without CPs but sim- 
ilar in other respects, none exhibit H2 emission when 
observed to similar detection limits. On the other 
hand Table 2 demonstrates that for a complete sub- 
set of Fabian's (1994) list of brightest cooling flow 
galaxies all CF galaxies with radio and optical emis- 
sion have been deteced in H2. 

In the detected galaxies the molecular gas is con- 
centrated within the central few kpc of the galaxy 
and on those scales is most likely co-spatial with the 
Ha-emitting gas. Since the close link between optical 
emission lines has already been established (Heckman 
et al. 1989), indicating that perhaps shocks may play 
an important role, and even correlated x-ray and op- 
tical filaments have been found in CPs (Sarazin et 
al. 1992), it is conceivable that also the H2 emis- 
sion is directly related to the CPs. Molecular line 
ratios are consistent with thermal excitation around 
T ~ 1500K. 

The theoretical interpretation for the connection 
between the cooling flow and the molecular hydro- 
gen is discussed in greater detail in JB97. These au- 
thors find that the molecular masses cooling through 
T~ 2000 K in the CF are a factor -100 too low to 
account for the H2 emission and suggest that continu- 
ous reheating of the gas by shocks or photoionization 
must take place. In this case only a small fraction of 
the total gas deposited by the CP is enough to ac- 
count for the emission. However, the high H2-to-Ha 
ratios we and JB97 find in the CP galaxies are dif- 



ficult to explain with current models. JB97 suggest 
that this requires H2 excitation by either slow shocks 
or suprathermal secondary electrons produced by X- 
ray photoionization of cold clouds. 

It is attractive to associate this excitation solely 
with the action of the CPs. However, all the galaxies 
detected in H2 have powerful radio jets that might 
be another ingredient underlying these observations. 
Por example, Seyfert galaxies too have narrow emis- 
sion lines with v ~ 400 km s^^ and sometimes strong 
molecular hydrogen emission. As shown by HST ob- 
servations (e.g. Capetti et al. 1996; Palcke et al. 1998), 
their emission line regions are often concentrated in 
distinct features (filaments and bow-shocks) indica- 
tive of a strong interaction between their radio jets 
and the interstellar medium. 

Ruiz, Rieke, & Shields (1997) claim that jet-powered 
photoionizing shocks are the most likely source of ex- 
citation for the H2 emission in Seyfert galaxies. Two 
Seyferts (in addition to the CF galaxy NGC 1275) in 
a sample of 28 showed ratios typical of CF galaxies: 
MCG 8-11-11 and Mrk 6. Ruiz et al. explain these 
observations by attributing the luminous H2 emission 
to excitation by jet-powered shocks as they interact 
with molecular material. The kpc-scale radio jets in 
Seyferts are often confined to the disk of the host 
galaxy, and so in some extreme cases may encounter 
a region rich in gas which provides an ideal work- 
ing surface for the jet to excite molecular hydrogen. 
It is therefore tempting to speculate that ellipticals, 
with their larger-scale jets and relative dearth of gas, 
require an external process such as a CP to trans- 
port large amounts of molecular material to a location 
where it can interact with the jet. This jet interaction 
scenario has been discussed in more detail already for 
PKS 1346-1-26 to explain the radio-optical alignment 
in this galaxy (McNamara et al. 1996). 

As an alternative one could imagine that the en- 
hanced H2 emission seen so far is a consequence of 
ongoing mergers in the cD galaxies rather than of the 
CPs. Further NIR spectroscopy of cD galaxies with 
and without CPs and high-resolution imaging with 
NICMOS will therefore be crucial to identify the ori- 
gin of the strong molecular hydrogen emission in these 
galaxies. 
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Fig. 1. — NIR spectra of galaxies with detected H2 lines in our sample plotted as Fx vs. A. The flux values are in 
detector units. Dotted lines give the estimated error derived from co-adding the individual spectra for each galaxy. 
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Table 1 
Emission Line Parameters 



Name 


z 


grating 


A 


flux 


EW 


FWHM 


Identification 






[l/mm] 


[Mm] 




[A] 


[A] 




PKS 0745-191 


0.1028 


300 


2.024 


(0.46) 


(31) 


(42) 


Ha 1-0 S(5) (A1.8353/im) 


[CF] 






2.068 


1.00 


55 


52 


Pao: 








2.086 


^0.06 


^3.5 


53 


H2 1-0 S(4) (A1.8914/im)" 








2.158 


0.50 


31 


55 


Ha 1-0 S(3) (A1.9570/im) 








2.166 


0.12 


8 


46 


Br7'' (or [SI VI]) 








2.241 


0.29 


19 


47 


Ha 1-0 S(2) (A2.0332/im) 


Abell 610 


0.0991 


300 


2.061 




<4 




Paa 


B2 0836+29 


0.0643 


75 


1.996 




<7 




Paa 


B2 0915+32 


0.062 


75 


1.991 




<5 




Paa 


PKS 1346+26 


0.06326 


75 


1.995 


1.00 


17 


(129^=) 


Paa 


(Abell 1795) 






2.014 


<0.18 


<3 


67 


Ha 1-0 S(4) (A1.8914/im)" 


[CF] 






2.081 


a56 


11 


72 


Ha 1-0 S(3) (A1.9570^im) 








2.162 


(0.24) 


(5) 


(52) 


Ha 1-0 S(2) (A2.0332/im)'* 








2.257 


0.66 


15 


81 


Ha 1-0 S(l) (A2.1213/im) 


PKS 2322-12 


0.0822 


300 


2.029 


1.00 


19 


39 


Paa 


(Abell 2597) 






2.048 


0.32 


7 


32 


Ha 1-0 S(4) (A1.8914^m) 


[CF] 






2.118 


2.02 


51 


50 


Ha 1-0 S(3) (A1.9570^m) 








2.200 


(0.63) 


(18) 


(46) 


Ha 1-0 S(2) (A2.0332/jm)'' 



Note. — (1) - Name ([CF] marks cooling flow galaxies), (2) - redshift, (3) grating used in observation, (4) - observed 
wavelength of line, (5) - flux relative to Paa, (6) - equivalent width in Angstrom, (7) - line width (FWHM) in Angstrom 
(not corrected for instrumental resolution), (8) - lino identification. Notes: "marginal detection, insufficiently corrected 
Br7 absorption in the caUbration star, '^broad fine width is most likely an artifact as indicated by the sudden increase of the 
noise in the spectrum, '^after subtraction of Br7 in the calibration star. 
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Table 2 

Brioiitest Cooling Flow Galaxies with Eadio and Hn Emission 



Name 


z 




M 


Ig P408MHZ 


L(Ha) 


Ha 1-0 






[10''^ erg/sec] 


[Mo/yr] 


[erg/sec/Hz] 


[10''° erg/sec] 




Abell 478 


0.0861 


241 


570 


31.4 


19" 


S(l)'' 


PKS 0745-191 


0.1028 


280 


702 


33.7 


490*" 


S(1)'*,S(2)-S(5)'= 


Hydra A 


0.0538 


35 


315 


34.2 


34== 


S(l)'* 


PKS 1346+26 


0.06326 


89 


478 


32.6 


184= 


S(l)-S(4)'= 


Cygnus A 


0.05654 


69 


187 


35.7 


2612= 


S(1)&S(3)^ 


PKS 2322-127 


0.0822 


64 


480 


33.3 


65*" 


S(2)-S(4)= 



Note. — Table of properties for galaxies from Fabian's (1994) list of brightest cooling flow galaxies, for which we have 
applied our original selection criterion {z > 0.0535, as well as radio and optical line emission) in hindsight. With the 
exception of NGC 1275 (redshift too low), this comprises the complete list of cooling flows with Ha detections published so 
far. The columns description is: (1) - name, (2) - redshift, (3) - 2-10 kcV luminosity, (4) - mass inflow rate, (5) - log of 
monochromatic radio power at 408 MHz (from NED database), (6) - Ha luminosity, (7) - Ha (1-0) lines detected so far. 
Additional notes: Columns 3 and 4 are taken from Fabian (1994), columns 6 and 7 were taken from: " White et al. (1994), 

Zirbel & Baum (1995), = Godon et al. (1994), JB97, = this paper, ^ Ward et al. (1991). We use Ho = 50 km sec"^ 
Mpc-^ 
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